M201 29 THE OPEN UNIVERSITY g 
Mathematics: A Second Level Course 
Linear Mathematics Unit 29 


Laplace Transforms 


J 


The Open University 


Mathematics: A Second Level Course 


Linear Mathematics Unit 29 
LAPLACE TRANSFORMS 


Prepared by the Linear Mathematics Course Team 


The Open University Press 


The Open University Press Walton Hall Milton Keynes MK7 6AA 


First published 1972, Reprinted 1976 
Copyright © 1972 The Open University 


All rights reserved. No part of this work may 
be reproduced in any form, by mimeograph 
or any other means, without permission in 
writing from the publishers. 


Designed by the Media Development Group of the Open University. 


Printed in Great Britain by 
Martin Cadbury 


SBN 335 01119 5 


This text forms part of a series of units that make up the correspondence 
element of an Open University Second Level course. The complete list of 
units in the course is given at the end of this text. 


For general availability of supporting material referred to in this text, 
please write to the Director of Marketing, The Open University, P.O. Box 
81, Walton Hall, Milton Keynes, MK7 6AT. 


Further information on Open University courses may be obtained from 
the Admissions Office, The Open University, P.O. Box 48, Walton Hall, 
Milton Keynes, MK7 6AB. 


1.2 


29.1 


29.1.0 
29.1.1 
29.1.2 
29.1.3 
29.1.4 
29.1.5 
29.1.6 


29.2 


29.2.0 
29.2.1 
29.2.2 
29.2.3 
29.2.4 
29.2.5 
29.2.6 


29.3 
29.3.1 
29.3.2 
29.3.3 

29.4 

29.5 


29.6 


Contents 


Set Books 
Conventions 
Introduction 


The Definition of the Mapping £ 


Introduction 

Piecewise Continuous Functions 
A New Type of Integral 
Functions of Exponential Order 
The Codomain of £ 

Can We Invert £? 

Summary of Section 29.1 


Some Techniques 


Introduction 

A Table of Laplace Transforms 

The Differentiation Formula 

The First Shifting Theorem 

The Second Shifting Theorem (Optional) 
Partial Fractions 

Summary of Section 29.2 


Applications of the Laplace Transform 


Differential Equations 
Electrical Engineering (Optional) 
Summary of Section 29.3 


Summary of the Unit 
Self-Assessment 


Appendices (Optional) 


1 The Convolution Theorem 
2 Green's Function and the Laplace Transform 
3 The Vibrating Spring: Impulse Functions 


Set Books 

D. L. Kreider, R. G. Kuller, D. R. Ostberg and F. W. Perkins, An Intro- 
duction to Linear Analysis (Addison-Wesley, 1966). 

E. D. Nering, Linear Algebra and Matrix Theory (John Wiley, 1970). 

Tt is essential to have these books; the course is based on them and will not 
make sense without them. 
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29.0 INTRODUCTION 


"This unit describes a method of solution for linear differential equations 
which is specifically designed to solve initial-value problems, and which 
is often more powerful for such problems than the more general methods 
we studied in Unit 9, Differential Equations II, Unit 11, Differential 
Equations ITI, and Unit 13, Systems of Differential Equations. In this unit 
we shall be more concerned with the rudiments of the theory than with 
actually solving difficult problems. 


The basic idea of the method is this. Suppose f is a bounded function in 
C[O, co). Its Laplace transform* is the function $ defined by 


$65) = [Pero di 
É 


o to 
where Í means lim [ » for all real values of s for which this limit 
o tolarge Jo 
exists. The mapping 


ifr 4 


with domain a suitable subspace of C[0, co) is a linear transformation, 
that is, a morphism with respect to linear combinations. The important 
thing here, however, is that under suitable conditions this mapping £ is also 
a morphism with respect to the unary operation of differentiation. It 
turns out for suitable functions f that the Laplace transform of f", the 
derived function of f, ist 


Ss [ero di  stL/X) —/0) 
If in addition /(0) — 0, i 

ijas [eo di = 508), 
where $ —£[/]. i 


Thus for a class of functions satisfying certain conditions the following 
diagram commutes: 


differentiate 


f 


multiply 


$ 


where sọ means the function 


sp 


by s is 


si sġ(s) 


Because of this morphism, the Laplace transform is of great value for 
solving differential equations, since it maps the operation of differentiation 
to the purely algebraic operation of multiplying by the function s ——— s. 
If f satisfies a differential equation, the Laplace transform will map this 
to an algebraic equation for $, and this algebraic equation is easier to 


* Marquis de Pierre Simon Laplace (1749-1827) was not only a great French mathematician 
and philosopher but also played a prominent part in French politics. He was, at one time, 
the Minister of the Interior under Napoleon, a post from which he was dismissed after a 
mere six weeks for " bringing the spirit of infinitesimals into administration” and elevated 
to a place in Senate. Although he is most celebrated for his work on gravitation it was in 
his analytic study of probability theory in the early I9th century that he developed the theory 
of Laplace transforms. 

f We write [f’] rather than (/7) to fit in with the notation of K. 
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solve than the original differential equation. For example, if f satisfies 
fi+feah 

and 
£0) =0 


where / is some given function, then @ the Laplace transform of f will 
satisfy 


sols) + é(5) = £C). 
This is an algebraic equation for $, whose solution is 
(5) = s + 1)" ' ED G)- 


To arrive at the solution of the original differential equation, therefore, we 
have only to calculate the inverse Laplace transform of @, i.e. the function 
whose Laplace transform is s ———+ (s + 1)^'£[A(s). There are well- 
established techniques for finding inverse Laplace transforms, and we shall 
study some of them in this unit. 
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29.4 THE DEFINITION OF THE MAPPING * 
29.1.0 Introduction 


We have seen above that the Laplace transform of a function fis given by. 
the formula 


f]: 5 —— [eo dt 
o 


but in order to define £ properly we should define its domain and codomain 
as well as its formula. In this section we look at these three elements of the 
definition. 


29.1.1 Piecewise Continuous Functions 


We want to be able to apply Laplace transform methods to as wide a 
class of functions as possible. As in Unit 22, Fourier Series, it turns out 
that we need not require f to be continuous; it is enough for it to be 
piecewise continuous. Indeed, one of the advantages of Laplace transform 
methods over the conventional method is that it requires no modification 
for piecewise continuous functions. To review the definition of piecewise 
continuity, 


READ Section 5-1, pages K177-8. 


Notes 


(i) lines 17-20, page K178 The restriction 1 € [0, co) describes the domains of 
the functions to which we apply the Laplace transform. All functions in the 
domain of £ have [0, co) as their domain. In the notation of Unit 22, the domain 
of £ is a subset of TCO, co). 

(i) lines 21-23, page K178 "piecewise continuous on every finite interval". 
For example, the square-wave function shown on page K177, which is periodic 
with period 2a, has an infinite number of discontinuities altogether; but it only 
has a finite number on any finite interval (of the f-axis) and therefore is piecewise 
continuous on [0, o). 


Example 


Consider the following electrical circuit. It contains an inductance L, a 


Vs 
A B 
x 
N 
==] 

V2 

AU N^ 

L R 


resistance R, and a switch that can connect either of two batteries of 
voltage V, and V;. The function 


Vir Vi) (te R) 


where V(t) = V, if the switch is connected to the ith battery at time t, 
describes the voltage between points A and B. If the switch is switched to 
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and fro, then V will be a piecewise continuous function. The differential 
equation describing the current i in the circuit is 


Li'(t) + Rit) = Vi) teR, 


and to solve this equation using the Laplace transform we would have to 
calculate 


[eo dt 


o 
which involves integrating a piecewise continuous function. 
Exercises 

i. Exercise 1, page K178. 

2. Exercise 2, page K178, omitting part (b). ' 


int 
(gm: iim “= 1) 


Solutions 
l. See page K739. 


2. (a) Piecewise continuous only if n = 1. This is because 


but 


does not exist for n > 1; so the limit of f(0 + A) ash 
tends to zero through positive values (see page K177, 
line — 6) exists only if n = 1. 


(c) Piecewise continuous on [0, co) (see page K178, lines 
21 to 23). There are an infinite number of discon- 
tinuities, but any finite interval [a, b] only contains a 
finite number of them (this number is less than b — a + 1). 


(d) Not piecewise continuous on [0, co), because there are 
an infinite number of discontinuities in the interval 
[0, 1]. (However, for any e > 0, the function is piecewise 
continuous on [e, 00).) 


(e) Piecewise continuous on [0, co), for the same reason as 
in (c). 


29.1.2 A New Type of Integral 


One of the new features of the integral defining the Laplace transform is 
its upper endpoint, which we have defined by 


o to 
f means lim Í FA 
o to large Jo 
provided the limit exists. An equivalent notation is 
to 
lim | f. 
toraJo 


If the limit does exist the integral Í fissaid to converge; if not, it is said to 
o 


diverge. 
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For example, 


EI z 
Í edi im | edt 
1 to^) 1 


= lim [—e"]9 
oe 


= lim (—e^'* + e`?) 


fortwo 


o 
Thus Í e~" dt converges to e~*. 
1 
On the other hand, 
o 1? to 
tdt= lim [5] 


= lim 402 —9) 


to^ 


which docs not exist. Ts | t dt diverges. 
3 


Example 

We calculate the Laplace transform of. 
ft!  (teR) 

Its transform £[f] maps s to 


EK) = fe" x 1 dt 


o 


" 
3 
LL 
| 

pa 
^ 
i 
x 
* 


Provided s > 0, we can make e^? as small as we please by making to 
large enough, and so 


lim (- D en =0. 
97m s 


That is 
o= 


if s > 0. It is typical of Laplace transforms that the integral converges only 
for s larger than some number, here 0. 


READ Section 5-2, page K179, as far as Equation (5-2) on page K180. 
Note 
line —4, page K179 f [cos ar] which stands for 
Ei ———» cos ar] 
is a function. 
Exercises 


I. Show that only one of the following integrals is convergent and evaluate 


it. 
(a) Í lg 
at 
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(b) L sin t dt 


o 
2. Find the Laplace transform (for s > 0) of 
fite—t (te R) 
(Hint: lim 1967? — 0 if 520. This will be proved in sub-section 
^^ 


t 
29.1.3.) 
3. Find £ [sin ars). 


Solutions 


1. Integral (a) is convergent: (b) is divergent. 


(2) Í ad= tim | lg 
zd maalit 


D 
S 
HE 
1 
ots 


1 ( 11,1 
-lim(-—]4z 
torw\ do 2 
zd 
-y 
o to 
(b) Í sin tde = tim Í sin t dt 
0 too JO 


= lim (1 — cos to). 
tono 
Since 1 — cos fọ oscillates perpetually between 0 and 2, 
the integral diverges. 


o 
2. w= [ e^t dt 
o 
1 
= lim [^(- LM 3) , by Section IILS.2 of TI 
toro s S7/1o 
t 1 1 
= lime*o( —-9—— = 
de ( s 3) M 
ls 
Haifs>0. 


o 
3. £ [sin at}(s) = f e`“ sin at dt 
o 


to 
= lim | e^*' sin at dt 
too 

to 


en 
= lim [ (—acosat —ssinat)| , 
to^ LS? + a? ) " 


by Section III.5.2 of TI, 


-sto 
= lim (5——35 (—a cos ato — s sin at, 
sig o o) 


fom 
+ 
Sa? 
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29.1.3 Functions of Exponential Order 


Tn order to define £ properly, we must specify its domain. We have already. 
indicated that this domain will be a set of functions that are piecewise 
continuous on [0, oo). We cannot, however, allow all such functions to 
belong to the domain of £, because we want the integral 


Í Safa dt 
0 


to converge for at least some value s, and there are functions (for example 
tr—— exp (?)) which increase so rapidly with t that the integral does 
not converge for any s. You will see in the next reading passage a simple 
condition of f which is sufficient (though not necessary) to ensure that the 
integral converges. 


READ from line 5 of page K180 to the end of Section 5-2 on page K182. 


Notes 


() line 1, page K181 L’H6pital’s rule does not form part of this course. It is 
a set of rules for finding limits such as 


MEC, 
Sims im Fey 


For our purposes the following will serve equally well. 


Theorem 


If n is a non-negative integer and a is a positive real number, then 


lim t"e-*' — 0. 
tam 


Proof. 


The Taylor series for the exponential function (Unit M100 14, Sequences and 
Limits II) is 


ER Uy, (ary 
Helter r ae 


(te R) 


1f at > 0, all the terms in the series are positive, so that 


(ary 
DU 


"> 


1t follows that 


men + 1)! 


ea 


Qi 0! 


qu 


re 


1 
x- 
t 


which can be made as close to zero as we please by making 1 large enough. 
We conclude that as / becomes large, t"e-*' must approach 0. 

(i) lines —12 to —9 page K181 The following diagram may help to make 
the comparison theorem plausible. The broken line shows the graph of f where 
it differs from that of |f}. 
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=t 


TOT 
0g 
Vy 
‘a 
Ta 


(iii) footnote, page K181 The comparison test for series is Theorem 1-6 on 
page K644. It was used in Unit 20, Euclidean Spaces I, sub-section 20.1.4. 

(iv) lines 2-3, page K182 “greatest lower bound ... all s «". You may 
prefer to think of so as the greatest lower bound of the set of real numbers s 
such that 


[ee 
o 


converges. For example, we saw in the previous sub-section (Equation 5-2) that if 
S(t) = cos at 


then the integral converges for s >0 but not for s <0. In this case the set for 
which the integral converges is *. Every negative number is a lower bound 
on this set, and so is zero; so so — 0 for -———+ cos at. 

A further discussion of greatest lower bounds can be found on page.K637 and in 
Unit M100 19, Relations. 

(v) line 7 page K182 "so may be — co". This means that in such cases the 
domain of C[/] is (— 2, co): the integral converges for all real s. 


Exercises 


Exercises 1, 2, 4, 5 on page K182. 


Solutions 
1 on page K182 
The Laplace transform (calculated in Exercise 2 of the pre- 
ceding sub-section) is 
1 
— z (s> 0) 
If s < 0 the integral diverges, and so the abscissa of convergence 
is 0. ? 
2 on page K182 


f[e"Ks) = lim j e^*te?! dt 


[E 


e(t ro 
= lim [ ] ifs#a 
Sto 


D 
8 


3 = 
foe d — S 
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Since the integral converges if s a and diverges for s < a the 
abscissa of convergence is a. 


4 on paye K182 
By the example in sub-section 29.1.2 
Ass 

EIG) = Zifs > 0. 
and the integral diverges if s <0. The abscissa of convergence 
is 0. 
5 on page K182 
In Exercise 3 of sub-section 29.1.2 we saw that 


£ [sin ats) =" 


zs 


The integral diverges if s < 0 and so the abscissa of convergence 
is 0. 


29.1.4 The Codomain of £ 


If £ is to be considered as a linear transformation, then its codomain must 
be a vector space. In this course we have dealt with many examples of 
vector spaces of functions. In most cases hitherto, every function in the 
vector space had the same domain. For example, every function in 
C[0, 1] has (0, 1] as its domain, and every function in C*(R) has R as its 
domain. 


The functions we obtain by applying the Laplace transform, however, can 
not all have the same domain. For example, we saw in the preceding sub- 
section that 


dió 1 
-stoot di = 
[^ e?! dt ae (57 a) 


o 


but that the integral diverges if s < a; so the Laplace transform of 


t-— get (t e (0, oo) 


se — (se(a, 00)), 


and the domain of the image function depends on a, that is on the function 
whose Laplace transform we take. This leads to some technical difficulties 
in defining a suitable vector space structure in the codomain of €; for up to 
now we have usually defined addition of functions when the functions 
added have the same domain. ‘ 


For the domain of £ we take £, the set of all piecewise continuous functions 
of exponential order. With the usual definitions of addition and scalar 
multiplication, & is a real vector space. For the codomain of £ we consider 
the set ¥ of all real-valued functions with intervals of the form (59, 00), 
[so 60) or (— co, co) as domain. To accommodate the fact that the mem- 
bers of F do not all have the same domain, we define addition in F as 
follows: if f, g € F, f+ g is the function whose domain is the intersection S 
of the domains of f and g, such that 


(f+ als) =f) + as), ses. 
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Having done this, we can see that F is not a vector space since there is 
not a unique zero element. For example if 


fx—1  (xell,o)) 
gxe—0 (xel, %)) 
h:xı—>0 (xe[0, o) 
then 
Jre=f 
and 
f*h-f. 


Since F is not a vector space, £: 6 ———» F cannot be a linear transfor- 
i i 
mation. Thus there exist f and g such that 


CU + g) s £L/] + Cle]. 
For example, let 
S(t) = cos at 
g(t) = —cos at 
By Equation (5-2), page K180 
Ei 


8f]: 5 — x (s e (0, o0) 
=s 
{[g):s > 12 (se (0, 20) 
so that 
CU] + tig]: s0 — (se (0, 0) 
but 
LU + eX = "47! x 0 dt =0 for all real s 
o 
so that 


t +g]:s 0 — (se(—00, «)). 


Thus £[/] + £(g] and £(/-- g] have different domains and are therefore 
different functions. 


READ from line 9, page K184 “ From this ...” to line 20, page K184.* 


Notes 


(i) lines 10 to 11 page K184 “CL f+ e] and £[/] +£[g] are identical for those 
values of s where both these functions are defined" means "f[/--g] and 
CLJ + Cle] map s to the same number for cach s in the domain of both these 
functions." 

Qi) lines 14, 15 page K184 “identical whenever they coincide on an interval 
of the form (a, «»)." This is'a similar device to the one used in Unit 22, Fourier 
Series, to define the Euclidean spaces SC[a, b]. We define a relation on the set F 
by the rule 


J~ g if there exists a real number a such that f(s) = g(s) for all s >a. 


It is not hard to show that this relation is an equivalence relation, and so we may 
say that f and g are equivalent if f — g, or as K says, the functions f and g are 
identical or identified. Note that different pairs of equivalent functions may 
coincide for different choices of a. For example, the functions 

$0——0 — (se(0, co) 
and 


sm0 — (se(—c, @)) 


* Note that we have omitted the early part of Section 5.3 of K in favour of the paragraph in 
this text preceding the READ instruction. 
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are equivalent (identical) since they both map all positive numbers to 0. Likewise 
the functions s ———» s and s ———,. [s| both with domain (—co, œ), are 
equivalent because every positive number s is mapped by both functions to the 
same image. Note carefully that we do not (indeed cannot) require that the 
images of the functions coincide on the entire intersection of their domains, 


S — isi 


at 


Gii) line-I7, page K184 “enforcing this identification”. This means that from 
now on K will write f= g (an abuse of notation) rather than f^ g to mean that 
the functions / and g are equivalent. 

(iv) footnote, page K184 This amounts to saying that K will write F to stand 
either for the space F* of equivalence classes or for F itself, according to 
context, In this text, however, we shall use ¥* for the space of equivalence classes. 


Exercises 
1. Determine the functions 
£[e'], £[e~'], t[e'] + £[e7*]. 


2. Which two of the following functions are equivalent under the relation 


~? 
ifs<2 
(s € (—09, œ)) 
22 
<ifs<3 
s 
g: s —— (s e (0, o2) 
n 
—dAidfísz3 
"nib 
lifs<4 
h: s —— (se (1, c2). 
1 
—Àifsz4 
s! s 


3. Prove that ~ is an equivalence relation. 
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Solutions 


1. From the exercises to subsection 29.1.3, we have 


£[e']: s —— ES (s& (1, ©) 


tle": s— — (se(-1, 05)) 
and so 
Ele] + £[e7]: s —— l pal (se (1, œ)) 
s-1 s+l 
since 


(1, œ) n (—1, œ) = (1, o»). 


2. We have 
as) = Ms) = — 3 for all s > 4 so that g ~ h 
= 
3. Reflexive 


* f ^ f" follows directly from the definition. 
Symmetric 


“If f~g, then g~f” also follows directly from the defini- 
tion. 


Transitive 

Now, if f ~ g, then for some ce R 
Sœ =s)  xe(oc) 

Similarly, if g ~ h, then for some de R 
g(x) =H) xeld o) 

Then if b is the larger of c and d we must have 
J) =g) =A) — xe(b co) 

and so f~ h. 


29.1.5 Can We Invert £? 


We have mentioned (in the Introduction) that in order to use the Laplace 
transform to solve differential equations we must be able to invert it; that 
is, given a function $ in F* we want a unique function f in & such that 
E*[/] = $. Whether this can always be done depends on whether £ is an 
isomorphism, i.e. whether £ is one-to-one and onto. These two points are 
discussed in the next reading passage. 


READ from line 21 of page K184 to the end of Section 5-3. 


Notes 


G) line —6, page K184 "differ only at their points of discontinuity”. Two 
functions f and g with the property mentioned are equivalent under the equiva- 
lence relation we used to define SC{a, 5] in Unit 22, Fourier Series, (page K332). 
We write this as fpg. For example, the functions 


fite——0 (rei, &)) 
and 
k ifreZ 
git 


0 if not Cei, 2) 
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are equivalent under this equivalence relation, so that Jpg. This is an equivalence 
relation on 6, and has nothing to do with the equivalence relation on F defined 
in sub-section 29.1.4. We will write &* for the set of equivalence classes under p. 
(Thus the mapping under discussion is really 


£*:8* —— $*) 
(ii) Lerch's theorem, page K185 This can be stated in terms of the two equiva- 
lence relations we have discussed: 

if f, g c6 and f[/] ~ fig], then fog. 


Gii) Last paragraph of the discussion on Theorem 5-2, page K185 The con- 
clusion is that £* is one-to-one but not onto; it is therefore not an isomorphism 
of &* to F*. If £-'[9] exists, it is unique (i.e. a unique equivalence class in £*,) 
but it does not exist for every $ in F*. Thus £-' is not a function on F*. 


Example. 


We have seen in Exercise 1 of sub-section 29.1.4 that 
1 


1 
Lie + e^]: s ———. 
$T 


(se (1, 09). 


s 
It follows that 


s 


or, more precisely, that 


Je erer 


£^ [54] is the equivalence class 
(: fe 6 and fp(t —— é + e7?)). 


Exercise 


Using the solutions to the exercises of sub-section 29.1.3, find the inverse 
images: 


(i) £-'[g] where 


1 a 
Satra Cte 


as) = 
and where a is a real positive number; 
(i) CHA], where 


I 
hls) = —— a<s<b 
Hs) s-a ( ) (se (a, «)) 
-0 (s>b) 
and where a and b are real numbers with a « b. (Don't jump to 
conclusions. Read the question carefully before trying it.) 


Solution 
(i) We know that i 

£fe") (s>a) 

and 

" a 

£[sin at] = mp (s > 0). 

Thus 
£fe" + sin at) = l + m 3 (sa) 

s-a s-a 
and 


17 


is the equivalence class in £* containing the function 
fpe e" +sinat (te [0, oo). 


This equivalence class consists of all real functions with 
domain [0, oo) which map only a finite set of points t in 
each finite sub-interval of [D, oc) to numbers other than f(t). 


(ii) Since A(s) = 0 for s > b, h is in the equivalence class in 7* 
which also contains the zero function. Since, by Lerch's 
theorem, £ is a one-to-one transformation from &* into F*, 
it follows that £^'[/i] is the zero equivalence class in &*. 
In other words, £7! [//] is the equivalence class of all functions 
with domain [0, co), which map only a finite set of points 
in each finite sub-interval of [0, oo) to numbers other than 
zero. The members of this equivalence class are the null 
functions (see page K332). i 


29.1.6 Summary of Section 29.1 


In this section we defined the terms 


Laplace transform (page K179) 

functions of exponential order (page K180) 

abscissa of convergence (page K182) 
Theorems 


1. (5-1, page K181) 


If f is a piecewise continuous function of exponential order, there exists a 
real number æ such that 


Í PN) di 
(J 


converges for all values of s > a. 


2. (page C11) 

If n is à non-negative integer and a is a positive real number, then 
lim t"e^*' = 0. 
t^ 


3. (5-2, page K185) (Lerch's Theorem) 


Let f and g be piecewise continuous functions of exponential order, and 
suppose there exists a real number s, such that 


CLs) = £leXs) 
for all s > sg. Then, with the possible exception of points of discontinuity, 
SO = 8 
for all t > 0. 
4. (5-3, page K185) 
If f is a function of exponential order, then 
lim £[/] — 0. 
se 
Techniques 


l; Find the abscissa of convergence of £[/] for suitable fe F. 
2. Find £[/] for suitable fe F. 


18 
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Notation 
ELIs) 
& 


F 
F* 


g 


(page K179) 
(page C13) 
(page C13) 
(page K184) 
(page C17) 
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29. SOME TECHNIQUES 
29.2.0 Introduction 


We have seen that calculating the Laplace transform of a given function f 
is a fairly straightforward procedure: one has to evaluate the integral 


pe 
| e^ f(t) dt. 
o 
Obtaining the inverse transform of a given function ġ is a more difficult 
matter; in principle, it amounts to solving an integral equation, that is, 
finding a function f such that 


4) = [e a (s€ (so, 0). : 
F 


In practice, however, various techniques are available which simplify this 
problem. In this way Laplace transforms become an efficient method for 
solving initial-value problems. In this section we discuss some of the tech- 
niques which are useful in inverting Laplace transforms. The situation is 
very similar to that which arises when we proceed from formulas for 
differentiation of functions to formulas for integration. We set up a table 
of derived functions, with as many functions as we care to list in the left- 
hand column, and the derived functions in the right-hand column. If the 
function we wish to integrate can be expressed as a linear combination of 
functions that appear in the right-hand column of the table (c.g. if it is a 
polynomial), then the integration is easy. If not, try the techniques of in- 
tegration by parts and by substitution, to knock the function into shape. 
It might work, or it might not. 


We are in a similar position with regard to £ and £^. If we wish to cal- 
culate £^'[f] for some f, the first thing to do is to try to express f as a 
linear combination of functions that we know to be images under £, If this 
fails, there are methods, analogous to integration by parts and by substi- 
tution, which extend the range of functions f which can be dealt with. 


29.2.1 A Table of Laplace Transforms 
The first step is to list some common Laplace transforms. 
READ Section 5-4 from page K186, to line 3 of page K187. 


Note 


Equations (5-7) to (5-10), page K186. We have derived all these in Section 
29.1. 


The table of transforms on pages K228-9 is important and will serve as a 
handy summary of the main properties of Laplace transforms. Ignore lines 
12 and 13 on page K229 unless you studied the delta function in the op- 
tional part of Unit 12, Linear Functionals and Duality. With this table, we 
will solve a simple differential equation by the Laplace transform. 
Example 
Find the solution of 

vy()—ay()me" — te [0, co) 


with the initial condition y(0) = 0. We will use the result quoted at the 
beginning of this unit, that if y(0) = 0, then 


£s) —s£D)G)- 


Its proof is given in the next sub-section. 
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Taking the Laplace transforms of both sides of the differential equation, 
we find (using the linearity of £) 


SLI) — ally Ils) = tfe7" yo) 
1 
O16) ~ ao 
(s— otis) e ——. 


Solving this algebraic equation we find 


1 
1 Baers} 


l a 
“a ( -a 
= tankot ( 
=; inh a Ks) 
by line 10 of the table on page K229. 
Thus 
Ties 
y) = z sinh at, 
by Lerch's theorem (since y must be continuous). 


Exercises 
l. Using the table on page K229, find: 


. a[2,5,.c 
() t IEEE 


(ii) ojh 2 3 ] 


s=] 8 5-3 


giving your answers as continuous functions. 
2. Use the Laplace transform and the table on page K229, to solve 
v(t) = ay() = e" 
y(0) = 0. 
Check your solution. 


Solutions 


: b fl m | axe ie | 
Lo eESS-« B] +00 ‘alte [5] 


ê 
cU Mop Doce. 


by the formula for £[r"]. 


1 2 3 
wy eef 2 3 
Gp t [atts] 


eph] pen] 


-e 4 2e? 3e 


by the formula for [e°] 
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2. The Laplace transform of the differential equation is 


tb] — ay] = t(r —— e") 


s£[yG) — af[y](s) 
since 

£D'YG) = stl) 
if y(0) =0. 
Solving this algebraic equation gives 

QUA) e 
and so 

y(t) = te” 

protest 

by the formula for £ nl with n = 2. 
Check y'(t) = (1 + athe"; 

so y(0) = 0 and y'(t) — ay(t) = e“ as required. 


29.2.2 The Differentiation Formula 


In the example in the preceding sub-section, we solved a first-order differ- 
ential equation with y(0) — 0 using the formula £[y'(s) = s£[yJ(s). We 
may, however, want to solve equations with y(0) #0, and equations of 
higher order than the first (for which additional initial conditions will be 
needed in order to give a unique solution). The next reading passage gives 
the generalization of the formula 

L'I) = sth) 


which we need to deal with such problems. 


READ from line 4 of page K187 to line —6 of page K189. 


Notes 


(i) line —7, page K187 "vanishes at its upper limit" Since f is of exponential 
order we have 


lim e-*'e f(to) « lim e7*'? Certo 
tone 


tone 
=0ifs>a 
and similarly 
lim e-* f(to) = — lim e-*? Ce"! = 0, 
fone dom 
Therefore 


lim e-*'? f(to) — 0 as stated. 


Qi) line —6, page K187. “jump discontinuity” was defined in Unit 22 (page 
K330). 

Gi) line 9, page K189 We shall discuss the method of partial fractions in 
sub-section 29.2.5. To follow the example, all you need to do is to recognize that 


(se R,s#0,5#1). 


(iV) line —10, page K189 "... of exponential order " What K seems to besaying 
here is that the study of constant-coefficient linear differential equations shows 
that al! their solutions are linear combinations of functions of the type listed on 
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page K180, line —5. This is certainly true as far as the general solution of a 
homogeneous constant-coefficient equation is concerned. The situation for 
nonhomogeneous constant-coefficient equations is a little more involved. It 
can be shown that if p(D) is any constant-coefficient linear differential operator 
and / is of exponential order, then every solution y of the differential equation 


P(D)y =h 
is of exponential order, but we shall not do this here. 
(V) line —8 page K189“... unique continuous solution." Since any solution of 


the differential equation must be at least n times differentiable, it is certainly 
continuous. Since the differential equation implies [y] — d, any such continuous 
solution must be a member of £-'[£], and then Lerch’s theorem shows that this 
continuous solution is unique. ` 


Exercises 


l. Solve 
»()—ay(t)-e* — te[0 co) 
with the initial condition y(0) — 1. 
2. Given that 
YO - y'() + yD =sin3t — te(0, co) 
with y'(0) = 4, y(0) = 5, find t[y]. 


Solutions 


1. Taking the Laplace transform and using Theorem 5-4 we 
obtain 


1 
£D) — af[y](s) = SEE 


1 
sts) — 1 — a [y](s) = en 
Solving for £[y](s) gives 


G- ott) e —— +1 


fp) = 


: [sinh ar] + t[e*'] 


fI. 
- tiz sinh at + e] 
a 
Thus 
d. " 
y(t) 2 — sinh at + e”. 
a 
2. The Laplace transform of the differential equation is 


3 
EDK) + Ks) + £D) = 55. 
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Theorem 5-4 gives 
LBK) — 5s — 4) + GEDAG) 7 5) + EDs) 


s.9 
and so 


1 3 
IO m re?) 


Notice how the manipulation of arbitrary constants has been 
replaced by the process of looking up inverse Laplace Trans- 
forms in a table. 


29.2.3 The First Shifting Theorem A 
READ Section 5-5, from the beginning on page K193 to line 10 of page K194. 


Notes 


(i) Equation (5-21), page K193 This is the most useful form of the first shifting 
theorem from the point of view of calculating images under £-'; to calculate 
images under f, use Equation (5-20). Equation (5-21) tells us what happens if 
we shift the graph to the left or to the right in the codomain of £. 


(i) line 2, page K194 The technique of completing the square is very useful in 
dealing with quadratic polynomials that have no real factors. (See, for example, 
page K134, line 4; or Section I11.2.3 of TI.) 


In general, the method is to write as? + bs + c, where b? < 4ac, in the form 
by V 4ac—b* 
dfe i) + “aa ] 
Exercises 
l. Using the first shifting theorem and the table on page K229, calculate 


G) Cle sin 24) 
Gi) Cle?) 


(i) ee- x sin ai]. 


2. Calculate 
" 1 E 
£-1)______|, i) £-'}=—__], 
j l +2)? + 7 © [ +454 d 
" 1 1 
t) i aj]. 
e E =I + jl WE [ae ES 3] 
Solutions 
1. (i) Since 
$ 2 
£[sin 27] = Fam 
(page K229), we have 
2 , 
£[e' si — 
[e' sin 21] Gone’ 
by Equation (5-20), 
Gi) Since 


2 
jn) -5. 
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we have 


2 
[e's] 2 — — 
d (sx 
(üi) Since 

s 
@+ay’ 


sin a| = 


we have 


aay t os s+2 
efe 3t si |= 
T ET 


to eur 1 1 
Pa 78) | eran 
w EFFE] EHE JE 


by Equation (5-21), 
= e^?! sin t (page K229) 


m te .—9 guo o5. 
a) zx] ii uc 
-2 
2eMpa|s5 
ý [5:4] 


= e^ ?'(cos t — 2 sin t) 

= 1 1 
ere | | 

GB) 4s— 1 + il 5 [i * j 


HEC [p EE 


= 4e' sin 4t 


9s? — 6s +2 


M fla] IER 


" eig 1 
9s? +1 


= seme 3] 


= gel" sin 3t 
(Further practice: Exercises 1, 2, 15, 31, 32 on pages K200-1.) 
29.2.4 The Second Shifting Theorem (Optional) 
Although this sub-section is optional we would encourage you to read it 
unless you are very short of time. 
READ from line 11, page K194 to line 2, page K197, and then from line —5, 
page K198 to the end of page K199. 
Notes 
(i) Equation (5-23), page K194 We can think of the unit step function as 
follows. If the domain of a function f includes points r <0, their images under f 


are irrelevant to the Laplace transform. We can imagine the integral in the defi- 
nition of Í to be over the whole of #, provided we first multiply f by us, because 


f ^r) dre 1 ” ky) di 
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for any function k: 4? ——> 2. In particular 


EX) -| “yn di -[ e7"f(t)us(t ) dt. 
o 


More generally: 


A(t) dr = | k(r)ust) dt. 


Thus the unit step function is a device for altering an end-point in an " infinite " 
integral. 

(ii) Equation (5-26), page K195 You may find the following form of Equation 
(5-26) easier to use: 


LL) = ez "£r —— f(r + 2) 


provided f(r) = 0 for t <a. It tells us what happens if we Shift: the graph to the 
left or to the right in the domain of £. 

Gii) Equation (5-27), page K196 This is the most useful for of the second 
shifting theorem from the point of view of calculating inverse images. 

(iv) Example 4, pue K196 The function / has a different formula in the three 
intervals [0, 1], (1, 2] and (2, 29). We want functions fi, f; , fs such that 


Sahthth, 


and the f, have formulas simple enough for their Laplace transforms to be written 
down directly. We achieve this by making 


fü) 9 fi) + Alt) + f) 
re (0,1) 1=f1)+0 +0 
re(1,2] 2-r—4) 00 
te (2, eo) 0 — f) + f) +A) 


To satisfy the equation for t € (0, 1] we put fi(r) =1, and falt) = f(t) = 0. 
Then we satisfy the next equation by keeping /,(r) = 1, fs(r) = 0 and by putting 
in the appropriate formula for /2(r), namely 1 — t. Finally, to satisfy the last 
equation, we put —(2 — 1) for fs(r), and keep fi(r) = 1, Alt -t. 

(v) line 13, page K199 "the sum of the geometric series..." The series 
converges only if s 2 0, and so the domain of £[/] is 2+. 


Exercises 


1. Use the second shifting theorem and the table on page K229 to cal- 
culate £[/] where 


(i) S() =u,(t) cos t 


tifüct«l 
di) f(-(2-tif1«t«2 1 
0iftz2 


(Hint: see Example 4, page K196) 


-— [sin tifü ct «zx 
Gi) Hye rris 


2. Calculate 


o ej 


l 


9 errn] 
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3. Calculate f[|sin:|] by using Theorem 5-9. (The function 
t —> [sin t| is the “ rectified sine wave” which is also discussed in 
Unit 22, Fourier Series.) 


fisinti 


Solutions 
1. (i) Eleni) cos (Ys) = e^ "P[cos(t + z))(s) (Equation (5-26)) 
= e^" t[—cos r)(s) 
= em 
ars 
Gi) SM =A +A + HO 
where 
Al) = tuolt) 
Salt) = (2 — 20u(0) 
Sst) = (—2 Qus). 


Hence 


ELA) = EUA] + £05] + £05] 


(page K229) 


LH) = z *e7t[ — 2(t + 1)](s) 
+e“ Ei- + (t+ 2) (8) 
(Equation (5-26)) 
-3 + e(z) + (5) (page K229) 


1-2e74e7 


= 


sey. 


Gi) SO =A + AO 


where 
fi) =sint 
Silt) = —u,(t) sin t. 
So 
£[f] = ffsin t] + Ef- u(i) sin (] 
= Cfsin 1] + e^ *t(— sin (t + z)] 
= f[sin t] + e^ "(sin t] 
Le; 


1 
8/16) = (04) o 


LM 29.2.4 


2. () Since 


i 
~1]_ | ae! K229), 
E F j e (page ) 
we have 
£7 = 
ec " ane 7? 
(Equation (5-27)). 
(i) Since 
t£ | + u =sint (page K229), 
s+] 
we have i 


" l eines 
£ lec -e'sint 


(Equation (5-21)) 
and so 


EY ME aid om 
t Geil = u;(Ne?‘ sin (t — 3) 


(Equation (5-27). 
3. Theorem 5-9 gives 
(1 —e7*) £(|sin t| Xs) -[ e^" sin t dt 
o 
> [= sin t — cos N 
~ +l " 
by Section III.5.2 (ii) of TI 
HE ud 
“os 41 
so that 


£f[sin :|]] = 


29.2.5 Partial Fractions 
Earlier in the unit (page K189) we evaluated an inverse Laplace transform 
by making use of the equation 
J 
s(s — 1) 
which expressed (in abbreviated form) the function on the left as a linear 


combination of two functions whose inverse transforms are tabulated. The 
table of transforms on page K229 lists only a limited number of cases, but 


(1) 


any function of the form A where p and g are polynomials such that the 


degree of p is less than that of q can be reduced to a linear combination 
Of the listed cases by the method of partial fractions. We consider, in 
place of s(s — 1), any polynomial q whose n roots are all real and distinct: 


qls) = G-a)Ys—-2a):G-2a). 


Then in place of the linear combination on the right-hand side of (1) we 
can consider the linear combination of the functions 


1 1 
s-a s-a ° s-a, 


(SER, S a...) Q 
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This leads us to consider whether all functions of the form 
(polynomial of degree <n — 1) 
Ee G) 
a(s) 
can be constructed in this way. In fact they can, as the following theorem 
shows. 
Theorem 


For a given polynomial function q, the set V of all functions of the form 
(3) is a vector space of dimension n, and if q is of degree n and has n distinct 
roots a,...,4,, then the set of n functions listed in (2) is a basis for V. 


Proof 
The first part of the theorem follows from the fact that the set of all poly- 
nomial functions of degree & n — 1 is a vector space of dimension n. For 


the second part, it is sufficient to prove that the n functions listed in (2) are 
linearly independent elements of V. Suppose that 


“i 


PES 
s-a s—a, 
is the zero function. We find on multiplying by 
(s = as — a2)... (sa) 
that 
a(s — a2)... (S — a,) 9 
+ (s — aja — a3)... (5— ap) 
Ut&G-a)e(-an)], (GeRssa,...,a) 
is the zero function on this domain. Although a, is not in the domain, we 
can take the limit as s approaches a,; this limit is 
(a, — az)... (a, — ap), 


and since a, ..., a, are all different we conclude that a, = 0. Similarly we 
can show that «, = 0, ... , &&, = 0, and hence that the n functions listed in (2) 
are linearly independent. By Theorem 3.3 on page N16, they therefore 
constitute a basis. 


Example 


$44 1 
Express ;—— — — as a sum of partial fractions. 
Xpress Ta 65 m p: 


We first factorize s? + 5s + 6, obtaining (s + 2)(s + 3). A basis for the 
space V of all functions of the form 


(polynomial of degree « 1) 


5*4 5546 
is { l " =} We therefore have, for some real «, and a, 
s+2 s43, 
s+4 [^ [2] 


SXST6 542 643 


544 —— a3) as 2) 
S3.5sR6 — (5-243) — 
Equating the coefficients of s and J in the numerators (which we may do 
since (s, 1) is a basis for the vector space of all polynomials of degree « 1), 
we obtain 


T=a, +a, 
and 
4 — 3a, + 22; 
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so that a, = 2, x, = — l, and the partial fraction expansion is 


s+4 2 1 


The method can be generalized to cases where the nth degree polynomial 
q does not have » distinct linear factors. As before, let V stand for the 
set of all functions of the form 


polynomial of degree < n — 1 
q 


where n is the degree of q. Then V has a basis consisting of the following 
elements: 


(a) for each distinct linear factor of the form s — a in q(s), there is an 
element ! 


$—a 
in the basis; 
(b) for each factor of the form (s — a)* in q(s), there are k elements 
Lei cds 
s—a (s-cay' ""'(s-ay 
in the basis; 


(c) for each factor of the form s? + as + b in g(s), with a? < 4b (so that 
s? + as + b has no real linear factors), there are two elements 


1 s 
S +as+b s? +as+b 
in the basis; 
(d) for each factor of the form (s? + as + b)* in q(s), with a? < 4b, there 
are 2k elements 
1 l 1 
s +as+b (sas by (s? as + bY 
s s s 
s +as+b (F +as+b) O? (s? + as - by 


in the basis. 


We shall not prove this result. 
Example 


To express in partial fractions, we use rules (a) and (c) above. 


y=] 
(s + Ds? + 1), 
This gives 
s-1 8 a as 
FIEF sel P41 FFl 
LG + 1) + @ + a3 s) + 1) 
j (s+ 1s? +1) 


Equating the coefficients of s?, s, 1 in the numerators, we obtain 


a, +a =0 
“+= l 


a ta, Ej 
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Thus 
a = —1,05; —0,a4 — 1. 
Example 
: s +s+] 
Find £^! |—— —7— — | by usi i i i 
in le Diet 5] by using partial fractions. Applying rules (a) 


and (b) above, we look for c, æ, , æ such that 
s +s+l [3 LA LA 
G*276*3 5:42 C+D! GD 
LAG + 2)G + 3) + a(s + 3) + als + 2)? 
(s + 2) (s + 3) E 


Equating coefficients of s?, s and 1 in the numerators, and solving the 
resulting system of three simultaneous equations for a,, x, «3, we obtain 


a, = —6, 4, =3,a,=7. 


The required inverse Laplace transform is therefore 


“| +s+] ]- "D TE IRSE: 
(5--2)(5-3) -— 542 (5-2) (543) 
= —6e7?! + 3te7?! + Te", 


The type of function that we can expand in partial fractions—the ratio 
of two polynomials—is called a rational function. Here we have only con- 
sidered the case where the polynomial in the numerator has degree less 
than the one in the denominator. For Laplace transforms this case is the 
only one that arises, because of Theorem 5-3 on page K185. 


Alternative Method 


In the case where the denominator g(s) of the rational function has 
degree n and has n distinct linear factors, the calculation of the coefficients 
0,, 05, ..., &, In 


p) _ aw 


can be perfomed without solving a system of n linear equations, by means 
of the following theorem. 


"Theorem (Optional) 


If 
a(s) = (s — ay)... (S — a) 
with all the a; distinct, and p(s) is of degree < n — 1, then in 


p) x om ; 
qs) Ei (s-a) 
each g; is given by 
me P(ai) 
$ (a; = aa; = a2) -+ (a; = a) 
omitting (a; — a;) 
(e.g. 
m pla) i 
(a = a3 = a3)... (a, = a) 
us plan) ). 
"(an = a), — a) «++ (8, = an1), 
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Proof 


[ect] A tam) 


JEL su 


20-40-02; 0:0. 
On the other hand 


p(s) 


(s — ajp(s) 


als) (s-a 


D 


G-a) 


omitting (s — aj) 
This function is continuous at a;, and so its limit is 


plai) 
(a= 4,)" 


omitting (a — a;) 


Exercises 

1. Write down partial fraction expansions, containing undetermined 
numbers a, 5, ..., for the following functions. Do not calculate 
ULL TELE 


T 1 
O > 


ET 1 
G) esd 
(i) s(s? +s +2) 


s(s +1) 
9 rl 


2. Calculate a;, &; , ..., for parts (i), (ii), (iv) of Exercise 1 and so obtain 


the appropriate partial fraction expansions. (If appropriate you can 
use the alternative method.) 


3. Use the answers to Exercise 2, and the formulas on page K229, to cal- 
culate 


@ t7 [4] 
o t" Lcx] : 


(i) e lesa 


Solutions 
b O S E dg (case (a)) 
O gp tey o) 
Gi) wep eee ets 
(case (a)) 
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v) = 1 EE esas 
SG? +542) s s$rst2 4542 
(case (c)) 
s(s + 1) [7 as a ys 
v = as 
© re r PH E Gap 
(case {d)) 
(i) We have 
az 
stl 


Bringing the right-hand side to a common denominator 
and equating powers of s and 1 on both sides, we obtain 


1 =a(s + 1) + a(s — 1) 


% +a = l and aq; —«; =0 
so that 

4% =4,0,=4. 
The expansion is therefore 


i 


stl 


In this case g(s) has distinct linear factors and so the 
alternative method can be used. Here 


Ps) =, 9(s)=(s — DG + 1); 
so 


and 


stl 


It is a good check to try a few special values of s; e.g. 
here s = 2 gives 

1 1 
2x1 2x3, 


2, 
$- 


1 zy [2 + a 
S(s+2? s “542° (542) 


(ii) Since 


we have 
1 = a(s +2)? + azs(s + 2) + azs. 


Equating coefficients of s?, s, 1 gives 


Oy +a =0 
day + 2a, a =0 
da, =I 
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whence a, = 4, 2; = — 4. 35 = — 4, and 
1 1 1 1 
sty 4s 4542) As+2) 
Check: s = 1 gives 
poro 0d I1 9-3-2 
374 4x3 2x3 36 


(iv) Since 
1 a a as 
z ait z 2 ub > 
s$ +s+2) s s +s+2 sjsct2 
we have 
1 =a (5? +s +2) + as + as? 
whence 
a+ x-0 ; 
ata, =0 
2a, =1 
so that 


1 1 1 


(SESH?) 2s As? +s+?) 


s 
~ 25? + s+ 2) 
Checks diu electi 1 
heck: $= 1 gives 7 = aka Iri 
s 1 1 
Lo Ese et a rete 
* 0 [x] [5*2] di 


(This one can also be obtained directly from page 
K229, as cosh t.) 


Ex ine 1 = = 
(i) € ox]: u pe^? 


(u) a p 


=f} 1 Giaa 


2s A(s+H* +3) A(s+ 2-3) 
=} — e7"? /$sin /it 
— ye^ "(cos 3t — 4/3 sin J30) 
Dod emin len A 


29.2.6 Summary of Section 29.2 


In this section we defined the terms 


unit step function (page K194) 

partial fraction expansion (page C30) 

rational function (page C31) 
Theorems 


1. (5-4, page K187) 


Let f be continuous on (0, co), and suppose that f" is piecewise continuous 
and of exponential order on [0, co). Then 


£Lf']= stu] - /(07) 
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where 

fQ*)-lim f(). 

1-0* 
More generally, if f, f’,...,f-" are continuous for all t > 0, and if f™ 
1s piecewise continuous and of exponential order on [0, œ), then 
{y"]= st] —sf(0*) — f'(0*) 

ELF] = EL] — EAO) — 72 (9*) — + — fon, 
2. (5-6, page K193) (first shifting theorem) EE 
If £[/] = $(s), then 

Ele] = (s — a). 
3. (page C29) 
For a given polynomial function q, the set V of all functions of the form ae 

(polynomial of degree < n — 1) 

q 


is a vector space of dimension n, and if q is of degree n and has n distinct 
roots a, ..., a, then the set of n functions 


1 1 
a, s-a, S a, 


E (SER, s #a,,..., a,) 


is a basis for V. 

4. (page C31) 

If 
g(s) = (s = a) + (s — a) 

with all the a; distinct, and p(s) is of degree < n — 1, then in 
p(s) 
as) (s= a) 

each g; is given by 


rn 
| (a; -a )(a; — a) 7 (a; 2) 
omitting (a; — a;) 
Techniques 
1. Given the Laplace transform of a function, find the function. te 
2. Use Laplace transforms to solve a constant-coefficient initial-value t 
problem. 

3. Given a rational function, find its partial fraction expansion. Hn 
Notation i 


ut) (page K194) 
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29.3 APPLICATIONS OF THE LAPLACE 
TRANSFORM 


29.3.4 Differential Equations 


We have already seen some simple examples of initial-value problems for 
differential equations solved using the Laplace transform. For an example 
of a more complicated equation, 


READ Section 5-6, page K202 to line — 10 on page K203. 


Notes 


(i) lines 8-11, page K203 This manipulation uses Theorem 5-8, which we have 
not asked you to read. The same result can be achieved without using the theorem, 
as follows: n 


t-i 3(s + 2) ati 3(s + 2) 
C++] M [s + 2)? + 9» 


-mC-i s 


= 3e7" csin 3r (line —3, page K229) 


(i) line —13, page K203 “The student who feels..." A further advantage 
of the Laplace transform method is that it is very suitable for dealing with cases 
where the nonhomogeneous term of the differential equation has different for- 
mulas in different intervals, like the function illustrated in Fig. 5-7 on page K196. 
The Laplace transform method is particularly useful for dealing with 
systems of differential equations. Consider, for example, the system 


2yy(0 + ya) + y.)  2y:() 
yy) + yal) +) + y) = 

with y,(0) = 1, y2(0) = 0. 
Since this system is not homogeneous, the methods we studied in Unit 13, 
Systems of Differential Equations, are insufficient to solve it. The Laplace 


transform method, on the other hand, yields the solution just as easily as 
for a single equation. The Laplace transforms of the equations are 


Qs + Dep) + G + DLs) = 2 + : 


(+ DD) +64 NDA) e 13 


Solving this pair of algebraic equations for £[y, (s) and £[y; (s), we obtain 


ob =- 1 (s+1-4) =e? 
1 Ll. T 29] 
£y) ai) ^23 zr 
In terms of partial fractions 
ut tod -—- 1-4 
Thus 
n-tE |; ie yite 4 2t 
stl] 
and 
yy aj L a) i.e. yz: t xe '=l=t 
StI s 
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Exercises 


l. Exercise 4, page K204. 
(Use the formula 


25s S 3 a 10 "ical: 
(s? + I)e— 2)?  s-2'(s—2y sz) 


2. Use Laplace transforms to solve the system 


y A Ys 


=% +y 


J; Mant + Pn 


In = Yn 


with y,(0) = ++ = y,-1(0) = 0, y,(0) = c. Here 2 is some real number; 
nis a positive integer, and y, ..., y, are the unknown functions. 


Solutions 


l. Taking the Laplace transform of the differential equation we 
obtain (using Theorem 5-4) 


SLDS) — des s — MCL) — Hy) + EN) 


so that 

=16 + 3s +p 2 s 
Z5- 2) -2 toa + 105-2) 
The inverse transforms of the various terms are 


—16 + 3s e 10 + 35 " 
mom" g7 I (Equation 5-22) 


£p) = 


£^! 


x B 
[oe 


Fa 10 42824 
-2 (G-2y s+) 


4 3 
— zz Sin t + =z cost 


25 25 
Therefore 

10 3 3 10 
ECT 291. 2r 4 10 eot 

x = e(- 25/*3 ze Be +75 te 

fi TE 
785 25 
= Pett inte S008 


2. Fori-1,2,...,n—], 
£D] = Alb + £i] 


(s — DEDS) = Chi). 
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The nth equation yields 


(s — ALDya](s) = c. 
It follows immediately by applying the previous equation 
with i — n — 1, then n — 2, etc. that 


— — 
Ga 


£pi) = 


and 


cin ie 

"(1 = ———- 

MOT 
This is an extremely neat method of obtaining the result. 
Compare it with the one used in sub-section 13.2.2 of Unit /3 
which involves numerous arbitrary constants. « 


29.3.3 Electrical Engineering (Optional) 


As we have seen, the great virtue of the Laplace transform is that it enables 
us to replace constant-coefficient differential equations by algebraic 
equations. This is particularly useful in electrical engineering, since it 
makes it possible to analyse networks containing inductors and capacitors 
by elementary methods instead of solving differential equations. For exam- 
ple, to analyse the network shown, the method described in Unit 4, 
Differential Equations I, and Unit 11, Differential Equations II (page 
K170), leads us to formulate the following differential equations. 


in(t) L iU, L 


voltage s 
source 
v(t) 


Kirchhoff’s first law 
d 
El =i- 120 


Kirchhoff's second law 


L a + ELO z wf) t> 0 (left-hand loop) ) 
di 1 | a 
L E + Rix()— 00-0 (0 (right-hand ioo 


We then have a third-order system of equations to solve. If the initial 
values of ij, i, and g are 0, we may solve the equations by taking the Lap- 
lace transform as follows. Putting 


Q = £[gl, 4 = tfi], 7; = Eli), V — £v], 
we get 


sQ(s) = (s) — h(s) 


sLI(s) +5 Q(s) = V(s) o 


sLh(s) + RIS) — Š Q(s)=0 
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Eliminating Q(s), we can simplify this to 


3G) o6) 14) = mal 
1 Q 
SLRS) + Rh) + — (5) = ()) =0 | 


The great virtue of this method is that we can go straight from the circuit 
to the system of equations (2) without having to bother about differential 
equations at all. All that is necessary is to apply Kirchhoff's laws to the 
circuit, treating an inductor just as if it were a resistor with resistance sL, 


and a capacitor just as if it were a resistor with resistance T We can now 
s 


solve the algebraic equations (3) for (say) /; in terms of V, obtaining 
sCV(s) 

€€——— — 

(LC + ISLC sRC+ 1) —1 

It is now possible to determine /; , and hence its inverse transform i, for 

any V= £[v], that is for any input voltage function v. But this is not 

necessarily the way an electrical engineer would actually proceed. Electrical 

engineers usually base their design criteria on the case where the input 

(here v) is, say, a sine function, so that V is proportional to 


DG)- 


a 
+a 
for some a. The Laplace transform / is then a rational function, and can be 
evaluated by the partial fraction methods we have described. (There are 


various tricks for speeding up the calculations, which could be rather 
tedious by the methods described in Section 29.2.) 


Exercises 


1. Without writing down any differential equations, obtain formulas 
giving J(s), the Laplace transform of the current i(t), in terms of R, C, 
and V(s) the Laplace transform of the voltage v(t) produced by the 
generator. Assume the initial current in the circuit and the initial 
charge on the capacitor to be zero. 


voltage 
generator 


yilt) 


2, Without writing down any differential equations obtain /;(s) in 
terms of R,, R2, L, C, and V(s). Assume the initial current in the 
network and the initial charge on the capacitor to be zero. 

i(t) Ri R2 

TAN NA^ 


MEG 


voltage V 
generator 
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Solutions 


1. Kirchhoff's second law applied to the Laplace transforms 
gives (see the discussion of Equation (2) above) 


V(s) = RKS) + = Ks). 


Thus 


2. Kirchhoff's first law tells us that the downward current in the 
inductor is i, (1) — i; (f). The second law then gives 


Vis) = Ryh(s) + LG) — 12) 
0 = —sL(A(s) — I4)) + Ra Is) + x h(s) 
Solving for /;(s) in terms of V(s), we obtain 


13) = sLV(s) 


—S— ENS, 
(sL + rofs +Ri+ x) -sL 
29.3.3 Summary of Section 29.3 

No new terms were defined in this section. 


Technique 


Solve a nonhomogeneous system of equations using Laplace transforms. 
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29.44 SUMMARY OF THE UNIT 


Laplace transforms are related to a linear transformation £* whose 
domain &* and codomain ¥* are vector spaces of equivalence classes of 
functions, defined as follows. 


&* consists of equivalence classes of piecewise continuous functions of 
exponential order, from [0, 90) to R, under the equivalence relation 


Fog if f(t) = gti) except at a finite number of points in every. 
finite sub-interval of (0, o9) 


("Piecewise continuous” means that the function has at most a finite 

number of discontinuities in every finite sub-interval of [0, co). “Of 

exponential order” means that there are constants C and a such that 
1/()| < Ce ((20)) 

F* consists of equivalence classes of continuous real functions with 

domains of the form (a, co). The equivalence relation here is: 


f ^ g if there is a constant c in the domain of both f and g, such 
that f(s) = g(s) for all s > c. 


The Laplace transform of f is defined by 


Of Ood gen 


and £* is then defined to map the equivalence class containing f to the 
one containing £[/]. The integral always converges for se(c, co) for 
some real number c, whenever fe 8. 


Lerch’s theorem states that £* is one-to-one. Its image space is a proper 
subspace of F*, since 


lim C[/Ks) =0 


for all fe & The application of the Laplace transform to differential 
equations follows from the fact that, whenever f e & and f has continuous 
derived functions up to the (n — 1)th derived function and a piecewise 
continuous nth derived function, then 


PDI Is) = SEKS) — $710) = 9779) - -- — 0709). q) 


This allows a differential equation in f to be converted into an algebraic 
equation in £[/]. In particular, if we restrict f to have 


LO) =f'0) = + =f") = 0, 
then for any polynomial P of degree n we have 
LIPIDY Ys) = POEL Is) 


giving a commutative diagram 


f AD (Df 
E 
£ i 
PSL) 
£[/] PO) = {[P(D)f] 
Thus, to find a solution y(r) to the differential equation 
P(Dy =f, 
with initial conditions 
»(0) = »'(0) = +++ = yO) = 0, 
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we proceed as follows: 
£[P(D)y] = £L] 
POLKS) = EKS) 


£| 
tye =O 
f 


For other initial conditions, the more general Equation (1) can beemployed. 


The problem thus resolves itself to that of finding £^ 1[/]. We do this where 
possible by means of a table of certain standard Laplace transforms (page 
K229), and the two shifting theorems: : 


LLS) = CLs — a) 
£[u.(D)g(t — ahis) = e^ *L[g](s) (optional material) 
where v, is the unit step function specified by 


u{t)=0 t«a 
=] t>a 


In particular, using the first shifting theorem, the table of standard Laplace 
transforms, and the technique of partial fractions, it is possible to find the 


; , Q 
inverse Laplace transform of any rational polynomial function 22 where 
degree of p < degree of g. 


Laplace transforms are especially useful in solving systems of linear 
differential equations and are used extensively in solving network problems 
in electrical engineering. 


Definitions 
Laplace transform (page K179) 
functions of exponential order (page K180) 
abscissa of convergence (page K182) 
unit step function (page K194) 
partial fraction expansion (page C30) 
rational function (page C31) 

Theorems 


1. (5-1, page K181) 
1f f is a piecewise continuous function of exponential order, there exists a 
real number « such that 


f en 'f(1) dt 
o 
converges for all values of s > a. 


2. (page C11) 
If n is a non-negative integer and a is a positive real number, then 


lim t"e^^' =0. 


two 


3. (5-2, page K185) (Lerch's Theorem) 
Let fand g be piecewise continuous functions of exponential order, and 
suppose there exists a real number sg such that 


ELIs) = tigis) 
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for all s > so. Then, with the possible exception of points of discontinuity, 
FO) =el) 

for all t2 0. 

4. (5-3, page K185) 

If f is a function of exponential order, then 
lim £[/] = 0. 


so 
5. (5-4, page K187) 

Let f be continuous on (0, œ), and suppose that f" is piecewise continuous 
and of exponential order on [0, co). Then 


£L] = st[/] — f(0*) 
where 
f(0*) = lim fð). 
moe 
More generally, if f, f", ..., 7?) are continuous for all t > 0, and if /? 
is piecewise continuous and of exponential order on [0, cc), then 


£[/] = ef] - sf(0*) — 70%) 


£L) = SELS) — =O) — 715 (0) = ++ — f 9095. 
6. (5-6, page K193) (first shifting theorem) 
If £[/] = ¢(s), then 
Eler] = $s — a). 
7. (page C29) 
For a given polynomial function g the set V of all functions of the form 
(polynomial of degree < n — 1) 
SS 
g 


is a vector space of dimension n, and if q is of degree n and has n distinct 
roots dj, ..., a, then the set of n functions 


1 
, Yen 
a 5—4; s—a, 


(se Rs #a,,..., ap) 


Li 
is a basis for V. 
8. (page C31) 
If 


g(s) = (s— a)... (8 — a) 
with all the a; distinct, and p(s) is of degree<n — 1, then in 
KOMT 
a(s) ii (s— a) 
each a, is given by 


Play) 


omitting (a, — a;) 


Techniques 


Find the abscissa of convergence of £[/] for suitable fe F. 

Find £[/] for suitable fe F. . 

Given the Laplace transform of a function, find the function. 

Use Laplace transforms to solve a constant-coefficient initial-value 
‘problem. 

5. Given a rational function, find its partial fraction expansion. 

6. Solve a nonhomogeneous system of equations using Laplace trans- 
forms. 


oe ee 
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Notation 
ELIs) 
& 


F 
g* 
e 

uk) 


(page K179) 
(page C13) 
(page C13) 
(page K184) 
(page C17) 
(page K194) 
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29.5 SELF-ASSESSMENT 
Self-assessment Test 


This Self-assessment Test is designed to help you test your understanding 
of the unit. It can also be used, together with the summary of the unit, for 
revision. The answers to these questions will be found on the next non- 


facing page. We Suggest that you complete the whole test before looking at 
the solutions. 


l. Which of the following functions are 
(a) piecewise continuous on their domains 
(b) of exponential order 
(c) in the domain of £? 


; if t>0) 
so—is if ica (t € [0, co) 


e' if t>n 

d Q0 if Daren (reto, o) 
eè if 128 

hit A if 6<r<a} Cel 0) 


2. Give the Laplace transform of 
t> e'sint-Fe^'cos2t (re [0, «o)). 
3. Express the following as partial fractions 


" 1 
® ri 
2543 
(s + 2)(s? — 4s + 1) 
645-5 +75 — 5* 
s(s — 3s + 2) 


(i) 


(iii) 


4. Give the inverse Laplace transform of 


. 1 

) s+ POSU (s e [0, 00)) 

i) s — — (re (0, co)) (optional part of unit) 
Ss 


5. If 
Y'(0 +x) = sin t, 


with y(0) = 3 and y'(0) = 4, find £[y](s) and hence the solution of the 
equation. 


6. Solve the following initial-value problem using Laplace transforms 
y" — »'() — 6y() 2 3 +1- 1 
with »(0) = —1, »'(0) = 6. 


7. Given the system 
1 ce nnne 
4 aM) = 


find the expressions for the Laplace transforms of y, and y,. 


45 


If the current in the circuit shown is i(f) and the charge on the capacitor 
at time 0 is 0, write down (without writing down any differential 
equations) ¢[i](s) in terms of s and the numbers A, C and F. 

(This question is based on an optional part of the unit.) 


constant 
voltage 


i 


i(t) 
> 
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Solutions to Self-assessment Test 


1. f (a) Yes (b) Yes ^ (c Yes 
g: (a) No (b) Yes (c) No 
h: (a) Yes (b No (c) No 
s ae. i stl 


+a 
s-1 F411 Gel s4 
3 (i) 1 A B ic 


a= + 5+ 
sS-1) s s sel 


_ As(s + 1) + B(s 1) + Cs? 
= s*(s 4-1) 


This gives the equations 


A+C=0 
A+B=0 
B=1 
hence 
A=-1,B=1,C=1 
and so 
i-e do di 
s(+1) s sS scl 
2543 A B Cs 


Gi) 


GDC asa) 542) oath! (s-41) 


_ As? — 4s + 1) + B(s + 2) + Cs(s + 2) 
j (s + 2s? — 4s + 1) 


This gives the equations 


A + C=0 
-4A+ B+2C=2 
A+2B =3 


hence 
rite, Beith. Gis 
A-c—$sB-i$Cois 
and so 


25 +3 wel 20 
(sa-2)5?—4s--1) 13(s42) © 1369? — 45 +1) 


5 
m 13(s? — 4s + 1) 


Gi) 6*5—5 4780-5 A B C 

x s$S(—3642 sS sS s 
N TE. 
"5-3 s+2 


_ A(s — 3G + 2) + Bs(s — 3X + 2) + Cs" (s — 3s + 2) + Ds(s + 2) + Es*(s — 3) 
» s*(s - 3/5 + 2) 
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This gives the equations 


—6A =) (6 
-A-6B = og 
A— B-6C =-1 


B- C+2D-3E= 7 
C+ D+ E=-1 


hence 
A=-1,B=0,C=0,D=4,E=-2 
and so 
64+5-s7 +79 —5* 4 9 
s(s— +2) s 55-3 5642) 
1 tt. 4 


V agen = tatr 
(from Solution 3(i)) 


Hence 
mal =-l+tt+e™ 
s*(s +1) 


(ii) cE] = u(t) 


s?£[yYG 


+l 


SEYNE) — 3s — 4 + f[ 


241 


1 
£[y](s) = i sles it 35+ P 
hence 
(t) = 4 (sin t — t cos t) + 3 cos t + 4 sin t 
= $ sin t + (3 — 4) cos t 
s*L[y](s) — yO) — y0) — stLy](s) + »(0) — ett) 


sfpl 1 
POP S 


s?£[yKGs) + s — 6 — sf[y](s) — 1 — 6tpy)G) -$ 


645—547) -s5 , 
EDs) = 96-3643 
NEM 
COMO = -3* 35—5- X63 
(from Solution 3(iii)) 
Hence y(t) = —40? + $e? — 2e7? 


3 
sth) - 2+ tb) = 5 


st(y2](s) — 1 4£DYG) = 


st 
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8. 


Solving for £[y, (s) and £[y2](s) gives 


rog d 
tb» <a E- e j] 
and 
3 3-3 
tbl) 2 *- c [-mee i 
£pi) = 1 
Ric 
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29.6 APPENDICES (Optional) 


These appendices are notes to parts of Chapter 5 of K, which deal with 
aspects of the Laplace transform that are not part of our syllabus, but 
which you may find interesting. In particular, Theorem 5-11, page K214, 
gives a characterization of the Green's function for a constant-coefficient 
linear differential operator. 


Appendix 1 The Convolution Theorem 


Notes to Section 5-7, pages K206—209. 
(i) Equation (5-33), page K206 The integrand here is the function 
gc fe- (Eeo) 
Different values of t give different values for the definite integral; thus varying ¢ 
produces a function, and it is this function whose Laplace transform is 
si dss). 


(ii) Equation (5-35), page K206 The function described in note (i) is here given 
a name and a notation. £ is a morphism from & with binary operation * to F 
with the binary operation of multiplication of functions 


(ha) Seg 
£ £ 
multiply 
CLA, fle) file] — CU'* 2] 


The convolution product has applications in algebra as well as analysis. 

Gii) line 4, page K207 Multiple integration is not part of this course, but it is 
analogous to integration in one variable: instead of integrating over intervals of 
A, one integrates over regions of #. For any particular region in #?, there are 
generally several ways of describing how to restrict the pair (£, 1) to lie within the 
region. In particular, one can restrict £ first, then describe the restriction on £, 
in terms of ¢; or vice versa. 


Appendix 2 Green's Function and the Laplace Transform 


Notes to Section 5-8, pages K212-216. 

(i) Equation (5-47), page K214 Using the formula 
£[D'yXs) = s'£[yG) — s*7*»() — +» — y^7X0) 

for k — 1,2, ..., n, we see that, if y obeys the conditions in (5-47), then 
£[pCD)»Ks) = pGYLYYG) — 1 


since the leading coefficient of the polynomial p is 1 by Equation (5-43). 
Thus the equation 


Ly=0 


(i.e. in our notation, p(D)y = 0), becomes, when the Laplace transform is applied 
to it, 


pisiy) — 1 =0 


ie. 
tp) zs 
ENT 
and so 
E EN 
wry =f [xs] (0. 


Gi) line 3, Theorem 5-11, page K214 The "auxiliary polynomial" of the 


operator p(D), is the polynomial p. 
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Appendix 3 The Vibrating Spring: Impulse Functions 
Notes to Section 5-9, pages K218—225. 


The first part of this section on the vibrating spring is very straightforward. 
The notes are all concerned with the concept of an impulse function. 


(D footnote, page K222 More generally, a force F(r) acting in a given direc- 
tion (but dependent on time) imparts an impulse 


1 
FQ) dr 

da 

to the body between times ro and f. 

This impulse is again equal to the change of momentum of the body due to the 

influence of the force over this period. 

When we give a sharp blow to an object, we succeed in imparting momentum to 

it. In other words, 


1 + 
FOd-| Fdo 
tty E 
where the blow occurred between a and a + 7. But the interval [a, a + 7] can 
be narrowed down to be as small as we please, as long as it contains the instant a 
at which the blow took place. As we saw in sub-section 12.3.2 of Unit 12, Linear 
Functionals and Duality, this implies that F cannot be a function as ordinarily 
defined (though it can be given an adequate meaning, as a linear functional). 
Qi) Equations (5-63), page K223, and (5-66), page K224 In the notation of 
Unit 12, sub-section 12.3.2, (7) and &(r — a) are So(¢) and 8,(r) respectively. A 
lot of the working can be short-circuited at this point, by noting that in sub-sec- 
tion 12.3.2 we defined 8,(1) in such a way that 


f fa fo 
whenever a € [b, c]. Thus, by definition, 


tlt — ay) = CBN = [ ec d 
A 
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